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computer-modelling  of  the  unstable  resonator  with  continuous  gain 


n n 


1 


I.  INTRODUCTION 

This  research  had  as  its  main  objective  a formulation  of  the 
equations  for  propagating  electromagnetic  waves  in  a laser  medium  in 
order  to  study  the  properties  and  behavior  of  the  laser  beam  in  gas 
laser  systems,  such  as  the  CW-EDL,  that  have  substantial  gain.  This 
study  of  these  beam  characteristics  and  factors  affecting  them  is  use- 
ful in  the  establishment  of  design  or  operational  modifications  that 
could  lead  to  improved  far-field  intensities  for  the  laser  output  or 
to  explanations  of  unexpected  beam  degradation  in  terms  of  laser  media 
or  laser  resonator  characteristics. 

In  order  to  investigate  the  properties  of  amplifying  media  of  gas 

lasers  relating  to  beam  quality,  the  application  of  a rigorously- 

derived  paraxial  approximation^  to  the  wave  equation  for  propagation 

through  an  amplifying  medium  was  necessary.  Laser  beam  quality  in  the 

device  is  a critical  parameter,  along  with  atmospheric  propagation,  in 

determining  the  far-field  intensity  pattern.  Factors  important  for  beam 

quality  are  diffraction  due  to  finite  beam  size,  beam  distortion  due  to 

saturation  of  the  medium  gain,  and  refraction  and  defocussing  due  to 

a difference  in  frequency  between  the  laser  frequency  and  the  center 

2 

of  the  gain  line.  While  Bridges  has  shown  the  effect  of  saturation  on 
propagation  in  an  amplifying  medium,  the  other  effects  have  not  been 
studied  in  depth.  After  derivation  of  the  equations,  they  are  solved 
numerically  on  a computer,  and  results  of  gas  lasers  with  various 
choices  of  the  important  parameters  are  plotted.  The  shapes  of  the 
curves  are  then  used  to  determine  beam  distortion  as  a function  of  the 


various  parameters. 
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As  the  result  of  studying  these  curves.  It  is  found  that  for  lasers 
with  media  of  the  general  characteristics  of  the  CW-EDL,  these  factors  of 
saturation  and  detuning  are  of  minor  importance  where  the  system  is 
treated  as  an  amplifier;  however,  further  consideration  indicates  that 
treatment  of  the  system  as  an  amplifier  is  not  adequate,  and  in  partic- 
ular, diffraction  may  be  of  significant  importance  in  the  case  of  an 
active  medium  in  a resonator.  On  the  other  hand,  detuning  is  very 
likely  to  be  of  minor  importance  in  the  resonator  case,  as  even  though 
the  system  is  disturbed  sufficiently  to  shift  its  operating  frequency 
from  line  center,  it  will  quickly  readjust  so  that  the  detuning  is  small. 

After  investigation  of  the  properties  of  an  amplifier,  the  equations 
were  then  set  up  for  a typical  CW-EDL  system  configured  as  an  oscillator 
using  a positive-branch,  confocal  unstable  resonator.  The  laser  system 
was  assumed  to  have  cylindrical  symmetry  in  the  gain  medium  for  the 


purposes  of  this  study. 
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II.  EQUATIONS  OF  PROPAGATION 

Standard  practice  in  theoretical  work  on  propagation  in  laser 

3 4 

amplifiers  and  in  the  study  of  modes  in  spherical  resonators  has 

depended  on  certain  assumptions  and  approximations  that  lead  to  an 
inconsistency.  Specifically,  it  is  assumed  that  the  electric  field  is 
plane  polarised  in  the  x direction.  According  to  the  exact  Maxwell's 
equations,  this  then  means  that  the  field  is  independent  of  x. 

However,  when  the  resulting  equations  are  solved  in  the  paraxial  ap- 
proximation, gaussian  solutions  are  found  in  the  transverse  direction, 
which  is  indeed  in  accordance  with  experiment,  in  spite  of  the  fact 
that  the  transverse  derivative  of  the  field  is  shown  to  be  zero  by  the 
exact  equations.  The  paraxial  approximation  has  been  analyzed  to 
resolve  the  inconsistency  and  the  resulting  equations  then  were  solved 
numerically  to  yield  quantitative  results. 

For  waves  propagating  in  the  z (axial)  direction  one  may  write  the 
field  as 

t -?T  + ae,  i*il“  <V  W <2-» 

where  the  subscript  T indicates  the  transverse  part  of  the  field  and 

a is  a unit  vector  in  the  z direction.  Let  the  total  gradient  be 
z 

7 " ’T  + lx  (2'2> 

where  represents  the  transverse  gradient.  One  now  takes  the  curl  of 

rf  a 

both  sides  of  the  Maxwell  equation  curl  t,  *>  iwv^H  (2.3) 

and  using  the  relation  curl  H = -iwe  k t,  (2.4) 

o 

one  obtains  curl  curl  (w/c)^k  t-  (2.5) 

where  w is  the  angular  frequency  and  k = I—  j + i 

\£o' 


where  e is  the 


dielectric  constant  and  o is  the  specific  conductivity.  Using  (2.1) 
and  (2.2)  in  (2.5)  the  transverse  and  longitudinal  components  of  the 
field  are  obtained: 

I - '°Tz  \ 2 - 3?T 

\ ft  ■ rr  ' 21k  ai1 

oZ 


2 - 

+ It 


-(f) 


ll  (VT*  V + lk7T-?T-^Fz=fe  <FZ 


(2.6) 


(2.7) 


Consider  a beam  of  width  w in  the  transverse  dimension.  It  has 

o 

2 ?TT  2 

associated  with  it  a diffraction  length  i = kw  = - — w where  X is 

° o X o 

the  wavelength.  These  characteristic  dimensions  may  then  be  used  to 
scale  (2.6)  and  (2.7) 


Let  x = w £ , y * w n,  z=£t 

o o 


Since  in  the  problem  to  be  considered  w <<  1,  let 

o 

w i 

f r — « — — 

l kw 


Using  (2.8)  and  (2.9)  in  (2.6)  and  (2.7)  one  obtains 

n fen  . T7  A c2  9F  V , \ 2-  ,3  32F 


(2.8) 


(2.10) 


(2.11) 


where  F (r  ,z)  -*■  F (p.O,  F (r  ,z)  -»  F (p.O 
T I T Z T £ 


(2.12) 


and  V e a |-p  + a ~ 

x * 3(  y 3n 


(2.13) 


with  A and  M unit  vectors  in  the  x and  y directions, 
x y 1 

The  dielectric  response,  k,  of  a medium  with  n atoms  per  unit  volume 


may  be  written  as 


,(><  - ok7!  (cgAo)J 


l + ft  + i 


L + (C8/,(jJ)  m^u»lF|2) 


(2.14) 


(2.15) 


where  ft  = (to-oo  )/y  . is  the  fractional  detuning  of  the  frequency  from 
ab  ab 

the  center  of  the  gain  line,  is  the  background  linear  dielectric  con- 

2 2 

stant,  g is  the  small  signal  gain,  and  I = |f|  /F  is  the  field  intensity 
in  units  of  the  saturation  intensity.  Choosing 


k2  = (oj/c)2  k. 


(2.16) 


equations  (2.10)  and  (2.11)  reduce  to 

2 2 -•  3 

V (fV  • F„  + f -r — *•  + IF  ) -fV  F -f  t— 

T T T 3?  4 T T 3; 


- 2 if  — = (fgim)  F 

T 


(2.17) 


f3  |-  (V  • f ) + ifV  • F -f2V2  F = [1  + (f 2gHm ) ] F 

3C  T x t x it,  6 c 


(2.18) 


As  in  physical  problems  of  interest  f gl  = g/k  = the  gain  in  a 
distance  k 3 is  always  small,  one  may  obtain  a consistent  solution  by 


expanding  the  field  in  powers  of  f.  Using  only  alternate  powers  one  has 
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for  the  transverse  and  longitudinal  field  components 


Ft  - Ft 


(0)  + f2  f (2)  + 


F,  + fFr(1)  + f3  Fr(3)  + . 


(2.19) 

(2.20) 


It  may  be  noted  that  (2.18)  contains  no  zero-order  terms.  Also,  since 

2 

rr.  is  a function  of  | F j , it  must  be  expanded  in  powers  of  f.  On  using 
(2.20)  with  (2.18)  to  first  order  in  f 

(0) 


F ^ = iV  • F 

C T T 


(2.21) 


from  which  it  is  seen  that  in  first  order  a small  longitudinal  compon- 
ent of  the  field  is  present. 

Using  (2.19)  with  (2.17)  the  lowest-order  terms  are 


V 2?  (0)  + 2i  — r-— 


1 T 


where  m = m 
o 


3 F 


(0) 


3*; 


-(8*>"0PT 


(0) 


*■(0)  i 2 


)=  ~-2- 
1 i + sr  + i 


(2.22) 

(2.23) 


(0) 


This  resolves  the  inconsistency  since  to  lowest  order  the  field  is  trans- 
verse and  may  depend  on  the  transverse  coordinate. 

Plane  polarized  solutions  of  the  form 


f (0) 


Ee 


iS 


T 

may  now  be  sought  where  E and  the  phase  S are  real.  Substituting 
in  (2.22)  and  equating  real  and  imaginary  parts,  one  has 


Vt2  E -(VtS)2  E-2  E = -(g  l)  (Re  mQ)  E 

2 (VT S • Vt)E  + VT2  SE  + 2 (||  ) = - (g  SO  (Im  mQ)  E 


(2.24) 

(2.24) 


(2.25) 

(2.26) 


taking  the  scalar  product  of  both  sides  of  (2.25)  and  (2.26)  with  E 
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III.  NUMERICAL  SOLUTION 

In  cylindrical  coordinates  equations  (2.27),  (2.28),  and  (2.33) 


become 


/ds\2  1 /s s\2  /3 S \ . „ ^1  3 / 3/T\  . 1 &T 

(~ — I + ~T  l ■TX ) + 2 5-7  ) = g i Re  m + r—  ( P 5— — I + ~ 5 — — 5- 

\3  pyl  n2\3(>/  0 m/T  3P  \ 3P  / pVl  3 $2 


13  / rJsK  1 3 lT  3 \ , 3 I _ » T - 

p 3p\Ip3pj  p2  3 <J>  I1  3 <J>  / 3c~_g£  mo 

dp  ^ 3 S 2 d<J>  _ 3 S 
d?  = fp  ’ p dC  “ 3 4> 

For  an  azimuthally  symmetric  amplifier,  these  equations  reduce  to 


dp  _ 3_S 

dc,  3 P 


i >-  M-U  |i  .-xu.. 

p 3 p \ 3 p / 3 £ o 


>2  + #) 


. „ . 1 3 

gi  Re  m + 

0 P/T  9 p 


I a/r\ 

(p  3 p j 


(3.1) 


(3.2) 


(3.3) 


where  p is  the  radical  distance  normalized  to  the  beam  size  wq,  C,  is 
the  axial  distance  normalized  to  the  diffraction  length  2. , I is  the 
beam  intensity  normalized  to  the  saturation  intensity.  S is  the  phase 
and  g is  the  small  signal  gain  per  unit  length  and  m^  * ^ J* - for 
a homogeneously  broadened  system,  fi  is  the  fractional  detuning  of  the 
operating  frequency  from  line  center. 

To  solve  these  equations,  we  introduce  the  beam  slope 


(3.4) 


and  use  it  to  eliminate  S in  Eqs.  (3.2)  and  (3.3) 


9 


p h (Ipv)  + H m " gi  1 Im  mo 


(3.5) 


3 , 2.  3v 

3T  (v  } + 2 H 


h (g*  ^ mo)+ [f"p  ~ptt  h[pfr)]  (3-6) 


with 


dl 

dc. 


3 I 3 1 

3 5 3p 


and 


dv  3 v 3_v 

dc"  3 £ V3p 


Eqs.  (3.5)  and  (3.6)  reduce  to 


~ = - g£  I Im  m - I-  - l|^ 

dC  op  3 p 


(3.7) 


dv 
d C 


1 3 p 

4 3p 


(3.8) 


where 


P=7g2+3-|  +?  + 2^Remc 


G - J-  (In  I) 


Eqs.  (3.4),  (3.7)  and  (3.8)  are  the  basic  equations.  They  can  be 
integrated  numerically  step  by  step  in  the  Cj-direction  along  the  path 
of  a ray. 

In  an  amplifying  system,  the  incident  beam  is  usually  given  at  the 
incident  plane,  e.g.  both  the  beam  intensity  and  the  beam  slope  are 
given  as  a function  of  p at  C “ 0.  Also  along  the  axis  of  the  amplifier, 
the  beam  slope  is  zero  and  the  beam  intensity  must  be  either  at  a 


J 
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maximum  or  at  a minimum,  e.g.  both  v « 0 and  C ■ 0 at  p = 0 for  all 
values  of  £.  With  these  initial  conditions,  Eqs.  (3.4),  (3.7)  and  (3.8) 
can  be  evaluated. 

The  numerical  procedure  used  is  to  first  divide  the  beam  at  the 

incident  plane  into  a number  of  rays  at  equal  distances  along  the  p- 

direction  and  integrate  p,  I,  and  v along  each  ray  step-by-step  in  the 

£ direction.  The  method  used  here  is  a Predictor-Corrector  Method. 

Assuming  that  p,  I and  v are  known  in  the  £ , and  C,  planes,  one  first 

n-i  n 

estimates  a tentative  value  of  Pn+1»  the  prediction,  in  the  £ ^ plane 


Vi  ■ + (If) 


where  is  the  integration  step  and 


M = v 
Mn  11 


from  Eq.  (3.4),  (simi- 


larly, for  the  other  two  quantities  I and  v) . One  then  calculates  a 


corrector 


where 


(tfiV 


v^+1  calculated  from  the  predicted  value.  If  the  dif- 


ference between  the  corrector  and  predictor  is  within  a small  pre- 
assigned tolerance,  one  stops  the  iteration  and  moves  to  the  next  plane; 
otherwise  the  calculation  is  repeated  using  the  corrector  as  the  pre- 
dictor and  recalculating  a corrector. 


This  method  is  used  throughout  the  integration  along  the  C direc- 
tion. For  the  first  integration  from  the  incident  plane  to  the  first 
plane,  the  extrapolated  values  Pn_j»  *n_i»  vn_^  were  used  to  calculate 


11 


the  first  predictors. 

In  evaluating  the  right  hand  sides  of  Eqs.  (3.7)  and  (3.8),  partial 
derivatives  of  various  quantities  with  respect  to  p are  required. 

These  are  obtained  by  fitting  a quadratic  curve  to  values  at  three 
nearest  points  to  the  point  where  the  derivative  is  needed.  The  deri- 
vative is  taken  to  be  the  tangent  to  this  curve. 

This  numerical  procedure  as  stated  can  be  applied  to  almost  any 
form  of  the  Incident  beam.  In  the  following,  consideration  is  restricted 
to  an  incident  collimated  Gaussian  beam.  Such  a beam  is  characterized 
by  two  parameters,  Iq  and  wq,  where  Iq  is  the  peak  intensity  in  units 
of  saturation  intensity  and  vp  is  the  beam  size  defined  to  be  that  at 

which  the  beam  intensity  falls  to  -7  of  I . Two  different  values  for 

e o 

w0  = 2.5  cm,  5.0  cm  and  three  for  1^  = 0.1,  0.25,  0.5  were  considered 
in  the  calculation.  In  a typical  calculation,  41  rays  are  taken,  the 
interval  between  successive  rays  at  the  incident  plane  is  taken  as 
0.05  wQ.  The  number  of  steps  in  integration  along  the  £ direction  are 
100  to  400  depending  on  the  length  of  the  amplifier  in  consideration. 

A sample  program  written  in  Fortran  IV  for  the  Univac  1108  system  is 
listed  in  the  appendix  1.  It  requires  approximately  20-30  seconds  for 
a single  calculation.  Figures  1-18  are  the  result  of  this  calculation. 
The  various  inputs  are 
Iq  = 0.1,  0.25,  0.5 
wq  = 2.5,  5.0  cm 
g * 0.01/cm 
Q » 


L 


0.0,  0.1,  0.5 
100,  200,  400  cm 
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where  L Is  the  length  of  the  amplifier.  In  these  figures,  the  amplified 
beams  at  the  output  plane  are  normalized  to  their  values  on  the  ampli- 
fier axis  to  make  the  distortion  more  evident.  These  normalization 
constants  for  various  inputs  are  listed  in  Table  I.  For  example,  for 
Iq  * 0.25,  wq  = 2.5  cm,  ft  = 0.5,  L * 400  cm,  the  amplified  intensity 
along  the  axis  at  the  output  plane  is  1.782  in  units  of  saturation 
intensity  or  an  amplification  of  7.13  times  of  the  incident  intensity. 

The  output  intensity  distribution  for  this  case  is  given  in  Figure  6 
as  the  outmost  curve.  In  Figures  1-18,  for  the  sake  of  comparison, 
the  input  gaussian  distributions  are  also  presented.  As  may  be  seen 
in  every  figure,  the  longer  the  amplifier,  the  larger  the  distortion 
away  from  the  gaussian  distribution.  The  amount  of  distortion  is  vir- 
tually independent  of  the  beam  size,  wq,  but  dependent  on  the  incident 
intensity,  Iq,  and  the  detuning,  ft.  The  distortion  is  large  for  large 
I and  small  ft.  All  these  features  may  be  understood  as  follows. 

Since  in  all  these  calculations,  the  length  of  the  amplifier  L is 
taken  to  be  much  smaller  than  the  diffraction  length  £,,  the  effect  due 
to  diffraction  and  refraction  is  neglibible;  i.e.,  the  beam  remains 
almost  parallel  to  the  axis.  This  was  also  confirmed  from  the  numer- 
ical calculation.  For  the  largest  deviation  in  all  the  cases  considered, 
a ray  starting  at  r = 2.0  wq  at  the  incident  plane  deflects  only  to 
r = 2.0005  wq  at  the  output  plane.  Thus,  if  one  neglects  the  effect 
due  to  diffraction  and  refraction,  Eqs.  (3.4),  (3.7),  (3.8)  would 
reduce  to  a single  equation 

dl  rKI 

d£  " 1 + ft*  + I 


(3.9) 
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IV.  UNSTABLE  RESONATOR  SOLUTIONS 

Although  the  solution  of  the  equations  for  a wave  propagating  In 

an  amplifying  medium  gives  some  Insight  into  the  effects  of  saturation, 

detuning,  and  diffraction  on  the  divergence  of  the  beam,  it  is  necessary 

to  model  the  amplifying  medium  in  a resonator  in  order  to  get  a more 

physically  correct  picture.  Diffraction  in  particular  is  likely  to 

assume  significance  in  the  case  of  the  resonator  due  to  repeated  passes 

and,  specifically  in  the  case  of  the  unstable  resonator,  the  presence 

of  the  edges  of  the  output  mirror  must  be  taken  into  account. 

The  most  applicable  work  done  in  the  past  on  unstable  resonators 

with  an  active  medium  was  that  of  Sziklas  and  Siegman^’^  and  that  of 

Rensch.^  Sziklas  and  Siegman  divided  the  laser  into  a number  of  axial 

segments  bounded  by  transverse  stations,  where  each  station  has  a fixed 

gain  and  phase  profile.  In  Ref.  5 they  carried  out  the  free-space 

propagation  between  stations  using  Hermite-Gaussian  functions  as  a 

basis  set  in  an  eigenfunction  expansion.  This  required  large  amounts 

of  computer  time  for  systems  with  large  fresnel  numbers,  which  are  of 

interest  in  high  power  lasers.  The  amount  of  computer  time  required 

4 

is  proportional  to  N , where  N is  the  fresnel  number.  In  their  later 

work.  Ref.  6,  they  increased  computational  efficiency  with  a discrete 

Fourier  transformation  carried  out  using  the  fast  Fourier  transform 

2 

algorithm.  In  this  approach,  computational  time  increases  as  N log^N, 
which  is  smaller  for  large  fresnel  numbers,  but  is  still  very  large 

j 

if  one  considers  fresnel  numbers  of  500  to  2000,  which  are  of  prime 
current  interest.  Rensch  used  a numerical  procedure  employing  a finite 
difference  method  to  solve  the  wave  equation.  The  resonator  was 

I 

I 
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divided  into  a series  of  axial  segments  with  each  segment  considered 

axially  uniform  in  saturated  gain  and  index  of  refraction.  Propagation 

through  each  segment  is  by  an  algorithm  that  satisfies  the  free-space 

wave  equation.  After  propagation  through  a segment,  the  amplitude  and 

phase  of  the  field  are  adjusted  according  to  the  average  gain  and  index 

of  refraction  for  that  segment.  Computational  time  for  this  approach 

2 L 

is,  for  a symmetric  resonator,  3N  -r—  , where  L is  the  length  of  the 

Az 

resonator  and  Az  is  the  length  of  the  axial  segment  over  which  the  gain 
and  index  is  considered  uniform.  Both  of  the  above  methods  have  a 
discrete  instead  of  continuous  medium  in  the  laser. 

The  previous  method  used  in  the  calculation  for  the  beam  distortion 
in  an  amplifier  system  is  now  extended  to  the  calculation  of  the  trans- 
verse mode,  beam  quality,  and  power-output  for  an  oscillator.  For  sim- 
plicity, the  oscillator  is  again  assumed  to  be  circular  with  axial  sym- 
metry, only  azlmuthally  symmetric  modes  are  studied. 

An  oscillator  in  general  consists  of  two  mirrors,  mirror  1 which 
will  be  called  the  output  or  the  front  mirror,  mirror  2,  the  back  mir- 
ror. These  two  mirrors  can  be  either  concave  or  convex  towards  to  the 
interior  of  the  cavity  with  different  radii  of  curvature  and  can  be  of 
different  sizes.  However,  for  most  high  power  lasers,  the  commonly 
adopted  system  is  usually  the  positive-branch  confocal  unstable  reson- 
ator. In  such  a system,  mirror  1 is  a small  convex  mirror  and  mirror 
2 a concave  mirror  with  much  larger  size  and  longer  focal  length. 

Figure  19  illustrates  such  a system.  Though  the  following  caiculational 
method  is  not  restricted  to  such  types  of  resonators,  the  emphasis 
throughout  is  on  this  type,  which  is  widely  used  in  large  laser  devices. 
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The  procedure  employed  for  calculating  the  transverse  mode  of  such 
resonators  is  to  introduce  two  planes,  plane  1 just  in  front  mirror  1 
and  plane  2 just  in  front  mirror  2 as  shown  in  Figure  19.  To  carry  out 
the  calculation,  an  initial  left-going  wave  with  its  intensity  and  slope 
is  assumed  at  plane  1.  It  propagates  to  plane  2 using  the  algorithm 
described  already  in  detail  in  the  previous  section,  since  propagating 
from  plane  1 to  plane  2 is  identical  to  propagating  from  the  input  plane 
to  the  output  plane  in  an  amplifier.  After  reaching  plane  2,  it  reflects 
back  from  mirror  2 to  plane  2 and  becomes  a right-going  wave.  This  wave 
is  propagated  to  plane  1,  again  using  the  same  alogithm.  After 
reaching  plane  1 and  mirror  1 the  beam  size  has  been  greatly  expanded 
and  is  now  much  larger  than  the  mirror  size;  therefore,  this  beam  splits 
into  two,  one  the  output  beam  passing  mirror  1,  the  other  a beam  reflec- 
ted back  to  plane  1 which  becomes  a left-going  wave.  This  calculation 
can  be  carried  out  indefinitely  until  the  wave  converges  to  a steady- 
state  stationary  distribution,  the  resulting  distribution  being  the 
single  transverse  mode  for  the  resonator  and  the  output  beam  passing 
mirror  1 is  the  output  of  the  laser. 

Presented  below  are  some  details  of  this  calculation: 

1.  In  propagating  the  beam  from  plane  1 to  plane  2 (and  from 
plane  2 back  to  plane  1),  the  gain  is  subject  to  saturation.  This  sat- 
uration gain,  unlike  that  in  an  amplifier,  is  determined  by  the  existing 
total  intensity,  l.e.,  the  sum  of  the  left-and  right-going  intensities 
at  this  point.  To  include  both  intensities  in  the  gain  distribution, 
the  intensities  of  the  on-going  wave  are  stored  in  memory.  These  values 
are  added  to  the  opposite  traveling  wave  intensity  on  the  next  pass. 
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Usually,  opposite  traveling  rays  do  not  cross  at  the  same  point,  in 
this  case,  a simple  interpolation  procedure  is  required  using  those 
stored  values  from  the  previous  pass  of  neighboring  points. 

2.  Reflecting  from  mirror  2:  Since  mirror  2 is  assumed  to  be 

larger  than  the  beam  size,  there  will  be  no  edge  effect.  The  field  is 
treated  as  though  it  were  physically  reflected  from  the  surface  of  the 
mirror.  Its  intensity  and  slope  are  adjusted  for  mirror  curvature, 
tilt,  and  reflectivity.  If  one  assumes  that  there  is  no  tilt  or  absorp- 
tion, to  a good  approximation,  there  will  be  no  change  in  intensity  and 
ray  position  as  the  distance  between  plane  2 and  mirror  2 is  extremely 
small.  The  only  change  is  in  the  ray  slope  with 

2p 

V = \r  — — I— 

refl  inc  R 

where  v^nc  and  v ^ are  the  slopes  of  the  incident  and  reflected  rays, 
respectively,  p is  the  ray  position,  R is  the  radius  of  curvature  of  the 
mirror  taken  to  be  positive  (negative)  for  a concave  (convex)  mirror. 

The  left-going  wave  now  becomes  a right-going  wave. 

3.  Reflecting  from  mirror  1:  In  addition  to  a change  in  ray 

slopes  upon  this  reflection,  mirror  1 is  smaller  than  the  size  of  the 
now  expanded  beam,  and  consequently  the  beam  upon  reflection  must  also 
be  truncated  according  to  the  shape  of  the  mirror.  At  this  point,  a 
severe  limitation  in  this  numerical  technique  is  apparent.  It  cannot 
treat  sharp  mirror  edges.  The  transverse  derivative  with  respect  to  p 
is  not  defined  at  the  mirror  edge  or  at  any  point  the  amplitude  is  dis- 
continuous. To  overcome  this  difficulty,  the  sharp  mirror  edge  is 
smoothed  out  by  introducing  an  effective  reflectivity  of  the  mirror 


! 
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which  has  a gaussian  profile  near  the  mirror  edge. 

R - e-(p-p0)2/"2  for  P £ p0 

and 

R = 1 for  p<PQ 

where  R is  the  reflectivity  of  the  mirror,  p^  is  the  mirror  radius  within 
which  the  reflectivity  is  uniform,  and  T is  the  truncated  distance.  For 
simplicity.  It  Is  further  assumed  that  the  distance  from  pQ  to  the 
physical  edge  of  the  mirror  is  equal  to  one  truncated  distance  T.  Fig- 
ure 20  shows  schematically  how  the  mirror  edge  is  treated.  Two  different 
values  T are  shown. 

4.  After  one  round  trip,  the  number  of  rays  reflected  from  mirror 
surface  1 is  greatly  reduced  as  many  of  them  will  pass  it  without  being 
interrupted  and  those  reflected  are  no  longer  separated  at  equal  dis- 
tance. In  order  to  regain  the  same  number  of  equally  separated  rays  as 
at  the  starting  point,  a simple  interpolation  procedure  is  used  based  on 
the  rays  reflected  from  mirror  1. 

5.  To  begin  the  calculation,  a left-going  wave  with  assumed  slopes 
and  intensities  is  launched  from  plane  1.  The  choice  of  these  initial 
values  will  make  a difference  in  the  amount  of  computer  time  needed  to 
yield  a stationary  solution.  A good  estimation  can  be  obtained  from  a 
simple  geometric  optics  calculation.  It  is  assumed  that  the  left-going 
wave  leaving  plane  1 is  a spherical  wave  whose  virtual  center  lies  at  a 

g 

point  behind  the  mirror  1.  The  intensity  is  assumed  to  be  a con- 
stant intensity  Iq  multiplied  by  the  effective  reflectivity  of  mirror  1. 
This  constant  intensity,  I can  be  estimated  from  the  gain  of  the  medium 
and  the  output  coupling  of  the  system  using  geometric  optics. 
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A specific  computer  calculation  is  made.  The  system  is  assumed  to 
be  confocal.  A complete  list  of  the  parameters  used  for  this  calculation 
is  listed  below: 


Radius  of  curvature  of  mirror  1 


Radius  of  curvature  of  mirror  2 


Length  of  the  cavity 


R^  = -535  cm 
R7  = 995  cm 
L = 230  cm 


Radius  of  mirror  1 


Radius  of  mirror  2 


Small  signal  gain 


a^  = 3.2  cm 
a^  = 6.0  cm 
g = 0.005/cm 


Frequency  shift  from  gain  center 


SI  = 0.0 


Truncated  distance 


T = 0.1  a. 


An  initial  spherical  wave  with  I = 0.50  I , where  I is  the  saturation 

oss 

intensity,  is  launched  from  plane  1.  Fifty-three  rays  at  equal  separa- 
tion, 0.025  a^,  are  taken.  The  number  of  steps  in  integration  along  a 
ray  from  one  plane  to  the  other  is  115.  After  seven  or  eight  return 
passes,  the  solution  converges.  The  total  computer  time  is  about  three 
minutes  on  a UN1VAC  1108  system.  Figure  21  shows  the  intensity  and 
phase  at  the  output  plane,  the  abscissa  scale  being  normalized  to  a^. 

The  vertical  line  at  1 represents  the  physical  edge  of  mirror  1.  The 
vertical  line  at  1.86  indicates  the  beam  size  at  the  output  plane  accord- 
ing to  geometric  optics  (the  output  coupling  for  this  system  based  on 
geometric  optics  is  71%).  Figure  21  shows  that  the  load  is  almost  uni- 
form on  mirror  1 and  the  output  intensity  increases  and  then  tapers  off 
very  quickly  near  the  geometric  edge  of  the  beam.  The  phase  is  also 
nearly  uniform  except  near  the  geometric  edge  of  the  beam.  Figure  22 


3hows  the  normalized  far-field  distribution  calculated  using  the 
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intensity  and  phase  of  the  output  beam  at  the  near-field.  The  far-field 

distribution  shown  in  Figure  22  is  very  nearly  identical  to  a diffraction 

limited  distribution  and  has  a first  side-lobe  peak  intensity  equal  to 

approximately  10%  of  the  main  lobe.  The  peak  intensity  at  0 = 0,  when 

2 

focussed  at  one  diffraction  length,  ka,  , is  found  to  be  1.21  I . The 

1 s 

dashed  curve  in  Figure  22  is  the  integrated  power  normalized  to  the 

total  power  (energy  in  the  "bucket").  As  shown,  more  than  80%  of  the 

power  is  within  an  angle  0 = and  more  than  90%  of  the  power  within 

/al 

0 = 1.5  x -x — . The  total  output  power  is  found  to  be  49.9  I watts 
-t  s 

1 2 
when  I is  measured  in  watts/cm  . 
s 

Figures  23  and  24  show  an  identical  calculation,  but  this  time 
with  ft  = 0.2.  The  result  of  this  calculation  is  almost  identical  to 
that  with  ft  = 0.0,  except  the  output  power  is  reduced  by  a small  per- 
centage. The  peak  intensity  at  far-field  is  now  1.14  I and  the  total 

s 

output  power  is  reduced  to  47.1  Ig  watts.  This  calculation  indicates 

that  the  effect  due  to  detuning  of  the  driving  frequency  from  the  gain 

center  is  not  important  as  expected. 

Figures  25  and  26  are  the  result  of  a calculation  with  the  small 

signal  gain,  g,  reduced  to  0.004/cm.  As  expected,  the  output  power  is 

greatly  reduced.  The  peak  intensity  at  far-field  is  now  0.63  I and 

s 

the  total  output  power  is  26.5  I watts.  The  intensity  and  phase  dis- 
tributions at  near-field  and  far-field  remain  unchanged. 

Lastly,  the  effect  due  to  mirror  edge  is  tested.  A different  trun- 
cated distance  T = 0.05  a^  is  used.  A smaller  T corresponds  to  a 

sharper  mirror  edge  and  is  closer  to  physical  reality.  Figures  27  and 

28  are  the  result  of  this  calculation.  It  is  seen  that  the  intensity 


and  phase  distributions  at  the  near-field  remain  almost  exactly  the  same 
as  before  except  near  the  geometric  edge  of  the  beam  where  they  fluctuate 
more.  The  far-field  pattern  is  also  almost  identical,  the  beam  being 
slightly  more  divergent.  The  total  output  power  remains  the  same,  but 
the  peak  intensity  is  down  to  1.07  I with  the  result  that  the  power  is 
down  about  10%  in  the  same  "bucket." 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

Laser  media  of  the  general  characteristics  of  the  CW-EDL  when  con- 
sidered as  an  amplifier  driven  at  a moderate  fraction  of  the  saturation 
intensity  show  little  beam  distortion  for  a gaussian  input  beam.  There 
are  only  small  saturation  effects,  and  very  little  effect  for  detuning 
of  the  driving  frequency  from  the  center  of  the  gain  line.  For  beams, 
such  as  those  considered  here,  which  have  diameters  of  many  times  the 
propagating  wavelength,  diffraction  due  to  the  finite  beam  size  is 
negligible. 

When  the  system  is  modelled  as  an  oscillator,  assuming  cylindrical 
geometry,  diffraction  effects  appear,  causing  a departure  in  phase  and 
intensity  from  the  uniform  phase  and  intensity  to  be  expected  on  solely 
geometric  considerations.  These  departures  from  the  uniform  plane  wave 
cause,  in  turn,  a change  in  the  far-field  intensity  pattern.  While  this 
azimuthally  symmetric  model  approximation  yields  some  valuable  insight 
into  system  behavior  and  is  of  great  value  in  establishing  model  and 
calculat ion  procedure  validity,  a model  in  rectangular  coordinates  pro- 
viding for  variation  in  the  gain  media  in  both  the  transverse  directions 
would  surely  be  more  accurate. 

One  observation  that  might  be  made  is  that  the  far-field  intensity  j 

pattern  might  be  Improved  by  increasing  the  system  feedback,  as  the 
medium  could  be  driven  harder.  This  could  be  accomplished  by  increasing 
the  diameter  of  the  output  mirror,  and  making  other  necessary  related 
changes.  Otherwise,  there  seems  no  advantage  in  changing  the  resonator 
geometry. 

I 

1 I 
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It  is  recommended  that  further  study  be  carried  out  on  this 
system  using  the  theories  and  procedures  developed  here.  Specifically, 
the  system  should  be  set  up  in  rectangular  coordinates  with  the  varia- 
tions of  the  media  in  the  transverse  directions  included  in  the  calcu- 
lations. The  effect  of  variation  in  feedback  on  the  far-field  intensity 
should  be  studied  on  the  computer,  as  it  is  relatively  complex,  inter- 
acting as  it  does  with  the  medium  saturation  and  the  output  aperture 
size. 


I 
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Input  Parameters  Peak  Intensity  at  Output  Plane 

in  Units  of  I. 

S 


w in  cm 
o 

— 
I ln  Ic 

O h 

r 

a 

L = 100  cm 

L = 200  cm  j 

L = 400  cm 

2.5 

0.1 

0.0 

0.237 

0.497 

1.436 

2.5 

0.1 

0.1 

0.235 

0.492 

1.419 

2.5 

0.1 

0.205 

0.392 

1.098 

2.5 

0.25 

0.519 

0.933 

2.115 

2.5 

0.25 

0.517 

0.927 

2.099 

2.5 

0.25 

0.5 

0.467 

0.798 

1.782 

2.5 

0.5 

0.906 

1.442 

2.783 

2.5 

0.5 

0.903 

1.435 

2.768 

2.5 

0.5 

0.5 

0.844 

1.301 

2.480 

5.0 

0.1 
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0.497 

1.436 

5.0 

0.1 

0.235 

0.492 

1.420 

5.0 

0.1 
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0.392 

5.0 
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0.519 

0.933 

2.115 

5.0 
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0.517 

0.927 

5.0 

0.25 

0.5 

0.468 

0.798 

1.788 

5.0 

0.5 

0.906 

1.441 

oo 

CNI 

5.0 

0.5 

0.903 

1.435 

2.771 

5.0 

0.5 

0.5 

1 

0.845 

1.303 

2.490 

TABLE  I Peak  intensity  of  the  output  beam  at  various  distances 
from  the  input  plane,  with  various  input  parameters. 
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FIGURE  CAPTIONS 


Fig.  1.  The  spreading  on  a gaussian  input  beam  at  various  distances 
from  the  input  plane.  The  output  peak  intensity  has  been 
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Fig.  20.  Mirror  edge  tapering,  T is  the  truncated  distance  as  defined 
in  the  text. 

Fig.  21.  The  near-field  intensity  and  phase  at  the  output  plane,  the 
solid  curve  is  for  the  intensity,  dashed  curve  is  for  the 
phase.  The  vertical  line  at  1 represents  the  physical  edge 
of  mirror  1.  The  vertical  line  at  1.86  represents  the  beam 
size  according  to  geometric  optics.  Input  parameters  are 
g = 0.005/cm,  ft  ” 0.0,  t=  0.1  aj. 

Fig.  22.  The  far-field  normalized  intensity  and  the  energy  in  the 
"bucket."  Input  parameters  are  g = 0.005/cm,  ft  = 0.0, 

T » 0.1  aj . 
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